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This paper shows a brief review on CDCC and the microscopic reaction theory as a 
fundamental theory of CDCC. The Kerman-McManus- Thaler theory for nucleon-nucleus 
scattering is extended to nucleus-nucleus scattering. New development of four-body CDCC 
is presented. An accurate method of treating inclusive reactions is presented as an extension 
of CDCC and the Glauber model. 

§1. Introduction 

The construction of microscopic reaction theory is one of the most important 
subjects in nuclear physics. It is a goal of the nuclear reaction theory. Furthermore, 
the construction is essential for many applications. Particularly for the scattering of 
unstable nuclei, there is no reliable phenomenological optical potential, since mea- 
surements of the elastic scattering are not easy. An important theoretical tool of 
analyzing inclusive reactions is the Glauber modelfQ The theoretical foundation of 
the model is shown in Ref. ED . The model is based on the eikonal and the adiabatic 
approximation. It is well known that the adiabatic approximation makes the removal 
cross section diverge when the Coulomb interaction is included. The Glauber model 
has thus been applied mainly for lighter targets in which the Coulomb interaction is 
negligible; see for example Refs. EJ) ,HJ) ,E|) JSJ ,d) JSD ,E!l and Refe-EHUm) for Coulomb 
corrections to the Glauber model. 

Meanwhile, the method of continuum discretized coupled channels (CDCC^P^'^ 
is an accurate method of treating exclusive reactions such as the elastic scattering 
and the elastic breakup reaction in which the target is not excited. The theoretical 
foundation of CDCC is shown in Refs. fT^|)J15pjl6p . Actually, CDCC has succeeded in 
reproducing data on the scattering of not only stable nuclei but also unstable nuclei; 
see for example Refs. ^^^^^^^^^^^^ and refer- 
ences therein. The dynamical eikonal approximation® is also an accurate method 
of treating exclusive reactions at intermediate and high incident energies where the 
eikonal approximation is reliable. The nucleon removal reaction is composed of the 
exclusive elastic-breakup component and the inclusive nucleon-stripping component. 
CDCC and the dynamical eikonal approximation can evaluate the elastic-breakup 
cross section, but not the stripping cross section. 

The experimental exploration of halo nuclei is moving from lighter nuclei such 
as He and C isotopes to relatively heavier nuclei such as Ne isotopes. Very recently, 
Takechi et al. measured the interaction cross section o~\ for the scattering of 28_32 Ne 
at 240 MeV/nucleon and found that a\ is quite large particularly for 31 NeP3 A 
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halo structure of Ne was reported with the experiment on the one- neutron removal 
reactionPD This is the heaviest halo nucleus in the present stage suggested experi- 
mentally and resides in the " Island-of- inversion" region. Determining the spin-parity 
of 31 Ne is essential to understand the nature of "Island of inversion". 

This paper shows a brief review on recent development of CDCC and the micro- 
scopic reaction theory that yields the foundation of CDCC. We present the micro- 
scopic reaction theory in Sec. [2] and new development of four-body CDCC in Sec. El 
We finally propose an accurate method of treating inclusive reactions as an extension 
of CDCC and the Glauber model in Sec. [H 

§2. Microscopic reaction theory 

In this section, we present a microscopic reaction theory for nucleus-nucleus 
scattering. This is an extension of the Kerman-McManus-Thaler formalism^ of 
the multiple scattering theory^ for nucleon-nucleus scattering to nucleus-nucleus 
scattering. In principle this reaction theory is applicable for many cases, but we 
consider the simple case in which the projectile breakup is weak, because it is not 
easy to perform fully-microscopic calculations including the projectile breakup. In 
the case, the theory is reduced to the double-folding model with the effective nucleon- 
nucleon (NN) interaction. The double-folding model is applied to nucleus-nucleus 
scattering at intermediate energies, particularly the scattering of Ne isotopes. The 
breakup effect is estimated by reducing the microscopic model to a three-body model 
and solving the three-body model with CDCC. This section is a brief review of 
Refs.EH,[35D,[36l). 

2.1. Model building 

Let us consider the scattering of projectile (P) on target (T). The most fun- 
damental equation for this case is the many-body Schrodinger equation with the 
realistic nucleon-nucleon (NN) interaction Vij. The multiple scattering theory^ £SJ 
for nucleon-nucleus scattering was extended to nucleus-nucleus scattering.® Accord- 
ing to the theory, the many-body Schrodinger equation is approximated into 

(K + h P + h T + Yl ^-£)^ (+) =0, (2-1) 
iePj'eT 

where E is an energy of the total system, K is a kinetic energy of the relative motion 
between P and T, and hp (/it) is an internal Hamiltonian of P (T). Here Tjj is the 
effective NN interaction in nuclear medium. The Brueckner g-matrix has commonly 
been used as r in many applications; see for example Refs. [37 |) .[38 |) .[39 |) .H0 |) . BT]) . 
Fi2"]) . l43p . l44p . l45p . The g-matrix interaction includes nuclear-medium effects, but not 
the effect of collective excitations induced by surface vibration and rotation of finite 
nucleus, since the interaction is evaluated in nuclear matter. The effect of collective 
excitations is small for intermediate energy scattering, as shown later. 

The Glauber model is based on the eikonal approximation for NN scattering 
and the eikonal and adiabatic approximations for nucleus-nucleus scattering. The 
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condition for the eikonal approximation to be good is 

\vnn(?nn)/Enn\ < 1, k NN a > 1 , (2-2) 

where Enn and fcjvjv represent a kinetic energy and a wave number, respectively, 
of the relative motion between colliding two nucleons, and a is a range of Vij. This 
condition is not well satisfied for the realistic NN potential that has a strong short- 
ranged repulsive core at small relative distance (rjvTv) between two nucleons; for 
example, ~ 2000 MeV at ttvtv = for AV18. 46> Actually, the eikonal approxima- 
tion is not good for NN scattering due to Uy, as shown in the left panel of Fig. [TJ 
To avoid this problem, a slowly-varying function such as the Gaussian form is used 
in the profile function of the Glauber modelP3 This procedure is justified in the 
following. 

In general, the (7-matrix has much milder r dependence than the bare NN poten- 
tial VijF* For example, the JLM (/-matrix® keeps this property. This means that 
the g-matrix is more suitable than v^j as an input of the Glauber model. Actually, 
the eikonal approximation is quite good for NN scattering due to the JLM ^-matrix, 
as shown in the right panel of Fig. [TJ The Glauber modeP is then more applicable 
to Eq. (|2Tp than the original many-body Schrodinger equation with Vij. In this case, 
the input of the Glauber model is not the profile function proposed in Ref. |4"T|) but 
the g-matrixPl 1 At higher incident energies where the Glauber model is used, the 
^-matrix is reduced to the i-matrix that has no medium effect. The t- matrix also 
has weak rjy jy dependenceP 1 This fact justifies the usage of the profile function with 
the Gaussian form. 
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Fig. 1. The on-shell NN scattering amplitude /nn((?) calculated with the bare NN potential AV18 
at the laboratory energy £?nn = 300 MeV in the left panel and with the JLM g-matrbP^ at the 
laboratory energy -Enn = 150 MeV in the right panel. The solid (dashed) and dotted (dash- 
dotted) lines show, respectively, the real and imaginary parts of /nn(<?) of the exact (eikonal) 
calculation. The left panel is taken from Ref. 

2.2. Application of double-folding model to reaction cross sections for Ne isotopes 

We analyze the scattering of Ne isotopes from a 12 C target at 240 MeV/nucleon. 
In the scattering the projectile breakup is weak, since the target is light and E is 
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large; see Sec. 12.31 for the point. The double-folding model becomes reliable in this 
situation. In the model, the potential U between P and T consists of the direct and 
exchange parts SHUHD The exchange part is non-local, but it can be localized with 
the local semi-classical approximation^ in which P is assumed to propagate as a 
plane wave with the local momentum within a short range of the NN interaction. 
The validity of this localization is shown in Ref. I50p . As the ^-matrix interaction we 
take the Melbourne interaction®'® that is constructed from the Bonn-B NN po- 
tential. 521 The projectile densities are constructed by either (I) by antisymmetrized 
molecular dynamics (AMD^S with the Gogny D1S interaction®'® or by (II) the 
deformed Woods-Saxon (DWS) modeP^'* with the deformation evaluated by AMD. 
Model I has no adjustable parameter, but the density is inaccurate in the asymp- 
totic region. Model II provides the density with the proper asymptotic form, but 
the model includes potential parameters. As the potential parameter set, we use the 
parameter set recently proposed by R. Wyss.® This set is intended to reproduce 
spectroscopic properties of high-spin states from light to heavy deformed nuclei, e.g., 
the quadrupole moment and the moment of inertia, and at the same time the root 
mean square (RMS) radius crucial for the present analysis. 




Fig. 2. (Color online) Reaction cross sections for the scattering of Ne isotopes on 12 C at 
240 MeV/nucleon. The experimental data for A — 28 — 32 are taken from Ref. I30[) . The data 
for 20 Ne is deduced from measured a\ at around 1 GeV/nucleoiP) with the Glauber model. ■' 
In the left panel, the solid (dotted) line represents the results of AMD (spherical Gogny-HF). 
A closed square is the result of AMD with the tail and breakup corrections. In the right panel, 
the dotted line represents the results of the DWS model, while the solid line corresponds to the 
results of AMD. 



Figure [2] represents the total reaction cross section ctr, for the scattering of Ne 
isotopes on a 12 C target at 240 MeV/nucleon. As shown in the left panel, the AMD 
calculations (solid line) succeed in reproducing the data,® whereas the spherical 
Gogny Hartree-Fock (HF) calculation (dotted line) undershoots the data; note that 
the spherical Gogny-HFB calculation yields the same result as the spherical Gogny- 
HF calculation within the thickness of line. The nuclei with A > 30 are unbound 
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in these spherical calculations. The enhancement from the dotted line to the solid 
line comes from the deformation of the ground state, since the deformation is a 
main difference between the two calculations. The deformation increases the or by 
at most 5%. The AMD results are consistent with all the data except 31 Ne. The 
underestimation of the AMD result for 31 Ne comes from the inaccuracy of the AMD 
density in its tail region. 

The tail problem is solved by the following resonating group method (RGM)P^ 1 
In principle the ground state ^( 31 Ne; 3/2J~) of 31 Ne can be expanded in terms of the 
ground and excited states <^( 30 Ne; J£) of 30 Ne. This means that the ground state of 
31 Ne is described by the 30 Ne+n cluster model with 30 Ne excitations. The cluster- 
model calculation can be done with RGM in which the ground and excited states of 
30 Ne are constructed by AMD. Here the wave function of 30 Ne includes many excited 
states with positive- and negative-parity below 10 MeV in excitation energy. This 
AMD+RGM calculation is quite time consuming, but it is done for 31 Ne. The tail 
correction to or is 35 mb that corresponds to 2.5% of <tr. The reaction cross section 
with the tail correction (a square symbol) well reproduces the experimental datcPS 
with no adjustable parameter. Thus 31 Ne is a halo nucleus with large deformation. 
The DWS modePO 1 well simulates the result of AMD+RGM for 31 Ne, as shown in 
the right panel of Fig. This may suggest that the DWS model is a handy way 
of simulating the AMD calculation with the tail correction. The difference between 
the AMD model and the DWS model for 28 ~ 32 Ne may show tail corrections to the 
AMD results. 

The same analysis is made for the scattering of Ne isotopes on a 28 Si target at 
38-60 MeV/nucleon. As shown in Fig. [3], the deformation effect is significant also for 
the lower incident energies, and consequently, the double-folding model yields better 
agreement with the experimental datap^ 1 where a normalization factor is multiplied 
in the theoretical results so as to reproduce measured or of 12 C+ 12 C scattering at 
38-60 MeV/nucleon. 
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Fig. 3. (Color online) Reaction cross sections for the scattering of Ne isotopes on a 28 Si target at 
38-60 MeV/nucleon. The experimental data are taken from Ref. l58|l . See Fig.[2]for the definition 
of lines. 
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2.3. Breakup, dynamical deformation and reorientation effects 

For a weakly bound system such as 31 Ne, the projectile breakup effect is not 
perfectly negligible. This effect is estimated by assuming the two-body model for 
the 30 Ne+n system and solving the three-body dynamics of the 30 Ne+n+ 12 C system 
with CDCC. Here the potential between 30 Ne and 12 C and that between n and 12 C 
are constructed with the double-folding model with the Melbourne g-matrix, and 
the potential between 30 Ne and n is made with the well-depth method; see Refs. I59|) 
for the potential parameters. The correction is 10 mb corresponding to 0.7% of ctr,. 

When the projectile is deformed in the intrinsic frame, the deformation enlarges 
the radius of the projectile density in the space-fixed frame and eventually enhances 
the reaction cross section. This static deformation effect has already been included 
in the double-folding model by making the angular momentum projection. Another 
effect is the dynamical deformation effect, i.e., the effect of rotational motion of 
deformed projectile during the scattering. This effect on ctr is small for intermediate- 
energy nucleus-nucleus scattering.® This has been confirmed with the adiabatic 
approximation to the rotational motion of projectile and the eikonal approximation 
to the relative motion between projectile and target! 3 -^ 1 In this subsection, the effect 
is estimated with no approximation. For this purpose, we consider the scattering of 
30 Ne from 12 C at 240 MeV/nucleon and do coupled-channel calculations between the 
+ ground state and the first 2 + state of 30 Ne. The projectile density is calculated by 
the DWS model with the deformation evaluated by AMD. The coupling potentials in 
the coupled-channel calculations are obtained by the so-called single-folding model. 
Namely, the nucleon- 12 C potential is first evaluated by folding the Melbourne-g- 
matrix interaction with the target density and the coupling potentials are obtained 
by folding the nucleon- 12 C potential with the projectile transition densities. 

In the single-channel calculation with no dynamical deformation effect, the resul- 
tant reaction cross section is 1469 mb. This result overestimates the corresponding 
result of the double-folding model by about 10 %, which is accurate enough for the 
present test. In the coupled-channel calculation, the resulting reaction cross section 
is 1468 mb. Thus the dynamical rotation effect on the reaction cross section is esti- 
mated as less than 0.1 %. The reason why the effect is small for intermediate-energy 
nucleus- nucleus scattering is shown in Ref . I34"|) . The integrated inelastic cross section 
to the first 2 + state is 2.9 mb. This is 0.2 % of <7r, indicating that o\ « or. 

The folding potential U is not spherical in general, when the spin of projectile is 
not zero. This reorientation effect is also tested by the coupled-channel calculation 
for the scattering of 31 Ne(3/2~) from 12 C at 240 MeV/nucleon, where the single- 
folding model is used. The resultant reaction cross section is 1512 mb, whereas the 
corresponding cross section is 1515 mb when the non-spherical part of U is switched 
off. The reorientation effect is 0.2 % and hence negligible for intermediate-energy 
nucleus-nucleus scattering. 

2.4. Summary on microscopic reaction theory 

We have constructed a microscopic reaction theory for nucleus-nucleus scat- 
tering, using the multiple scattering theory. This is an extension of the Kerman- 



Recent development of CDCC 



7 



McManus-Thaler theory for nucleon-nucleus scattering. In the theory, nucleus- 
nucleus scattering is described by a multiple scattering series due to the ^-matrix NN 
interaction instead of the realistic one. The former has much milder r-dependence 
than the latter and therefore the Glauber model is applicable for this theory with 
high accuracy. 

In the nucleus-nucleus scattering from lighter targets at intermediate energies, 
the breakup effect is small and hence the double-folding model becomes reliable. In 
this situation, we can use the double-folding model in which projectile and target 
densities are constructed with fully microscopic structure theories such as AMD, 
HF and HFB. These fully microscopic theories have been applied to measured reac- 
tion cross sections for Ne isotopes. In the Island-of-inversion region, the nuclei are 
strongly deformed. In particular, 31 Ne is a halo nucleus with strong deformation. 
The dynamical deformation effect and the reorientation effect are also found to be 
small, so that the interaction cross section is identical with the reaction cross section 
with high accuracy. 

§3. Four-body CDCC 

In this section we present a method of treating four-body scattering in which the 
projectile breaks up into three constituents. This method is called four-body CDCC. 
In this method, the three-body continuum states of projectile are discretized by di- 
agonalizing the internal Hamiltonian of projectile with the Gaussian basis functions. 
This discretization is called the pseudo-state method. Recently, we proposed a novel 
method of calculating continuous breakup cross sections with the complex-scaling 
methodP-ES This section is a brief review of Refs. 122]) .1231) .1281) . 

3.1. Formulation 

We consider the scattering of 6 He as a typical example of four-body breakup 
reactions. Here we focus our discussion on the differential breakup cross section. The 
scattering of 6 He on a target A is described by the four-body Schrodinger equation 

[H-E tot }\^} = (3-1) 

with the outgoing boundary condition, where the total energy E to t satisfies E to t = 
+ £q for the corresponding incident energy J5-? M in the center-of-mass system 
and the ground-state energy eo of 6 He. The total Hamiltonian H is defined by 

H = K R + hp + U nA + U nA + U aA + y a c A oul (3-2) 

with 

hp = K y + K r + V nn + V na + V na , (3-3) 

where hp is the internal Hamiltonian of 6 He. The relative coordinate between 6 He 
and A is denoted by R and the internal coordinates of 6 He are by a set of Jacobi 
coordinates, £ = (y,r). Momenta conjugate to R and (y,r) are represented by P 
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and (p, k), respectively. The kinetic energy operator associated with R (£) is repre- 
sented by Kr (K^), V xx i is a nuclear plus Coulomb interaction between x and x', and 
U x a and V^4 Ul are nuclear and Coulomb potentials between x and A, respectively. 

In CDCC with the pseudostate discretization method, the scattering is assumed 
to take place in a modelspaceP^'EU'EZI'ED 

7 

where ^ 7 is the 7-th eigenstate obtained by diagonalizing hp with L 2 -type basis 
functions. The four-body Schrodinger equation is then solved in the modelspace: 

V[H-E tot ]V\^ cc )=0. (3-5) 

The modelspace assumption has already been justified by the fact that calculated 
elastic and breakup cross sections converge with respect to extending the mod- 
el sp aceJ2DJZDi2aJ23» 

The exact T-matrix element to a breakup state with (p, k) can be described by 
T £ (p,k,P) = (4-\p,k) x i-\P)\U -Vf^\vM), (3-6) 

where V^° ul is the Coulomb interaction between 6 He and A. The final-state wave 

functions, \p,k)) and \ X i with the incoming boundary condition, are 

defined by 

Tr + ^ C h° c u1 - (£tot - e)] \xi~\P)) = 0, (3-7) 
[h P -s]\4-\p,k))=0, (3-8) 

where E tot — e = (HP) 2 / (2/xr) and e = {hp) 2 / (2p y ) + (hk) 2 / (2p r ) for reduced masses 
Pr and /i£ of coordinates R and respectively. Inserting the approximate com- 
plete set Eq. (|3-4p into Eq. (j3-6|) . one can get the T-matrix element with high 
accuracypS'ED'EM 

T E (p,k,P) « ^(4-)(p,fc)|^ 7 )r 7 (3-9) 

77^0 

with the CDCC T-matrix element 

T 7 = (<P, X (-\P,)\U - V 6 c n °: l \^ cc ) (3-10) 

to the 7-th discrete breakup state <P 7 with eigenenergy e 7 . Here Eq. ()3-9p is derived 
by replacing P by P 7 in xi \p)- The T 7 are obtainable with CDCC, but it is quite 
hard to calculate the smoothing factor (i/jg \p, fc)|# 7 ) directly with either numerical 
integration^ 1 or the complex-scaling methodP^ 1 Hence, we propose a new way of 
obtaining the differential cross section with respect to e without calculating the 
smoothing factor. 
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Using Eq. ()3-9p . one can rewrite the differential cross section into 
d 2 a 



j d P 'dk'5(e-e')\T £ ,{p',k',P')\ 2 « hz{e,Q P ) (3-11) 



dedQp 

with the generalized response function 

K(e,f2 P ) = Im\ ^ r*(# 7 |G^|# 7 /}Ty , (3-12) 

where G^ - -* = lim (e — hp — in) -1 . There is no smoothing factor in Eq. (13- 121) . as 

expected. The propagator operates only on spatially damping functions 

This makes the calculation of (^^G^ |^y) feasible. 

In order to calculate (# 7 |Gr(~) |^ 7 '), we use the complex-scaling method in which 
the scaling transformation operator C (9) and its inverse are defined by 

(r,y\C(9)\f) = e™f(re id ,ye w ), (3-13) 
{f\C-\e)\r,y) = {e~ m f{re- w ,ye- w )Y. (3-14) 

Using the operators, one can get 

(# 7 |(?H|# 7 ,) = ($ 1 \C- 1 (9)G < f ) C{e)\^), (3-15) 

where 

G^ = lim -s . (3-16) 

v^+o e - h p -irj 

with hp = C(9)h v C~ l (9). When -vr < 8 < 0, the scaled propagator {tlG^ \g) is 
a damping function of £ and note that is negative since G (_) has the incoming 
boundary condition. The scaled propagator can be expanded with L 2 -type basis 
functions with high accuracy: 



G (-) kE MH|1, ,347) 



where <f>f is the i-th eigenstate of hp in a modelspace spanned by L 2 -type basis 
functions, (4>i\h P \<f>f,) = efSu>. Inserting Eq. (|3Trl) into Eq. ^U\\ through Eq. 
(|3-15p leads to a useful form of 

d 2 a 1 ^ Tfff 



1 rpyrpy 

-ImV^, (3-18) 
>• , vr 4^e-ef 

with 

2? = ]T<^|C<(fl)|<2y>T y , If = ^^(^IC- 1 ^)!^). (3-19) 

7' 7 

This method does not require to calculate the exact three-body continuum states 
^ ( f\k,p). The convergence of this method is shown in Ref. 
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Fig. 4. Breakup cross sections for (a) 6 He+ 12 C scattering at 240 MeV/ A and (b) 6 He+ 208 Pb scat- 
tering at 240 MeV/A. The solid lines are results of full-fledged the four-body CDCC calculations. 
The dot-dashed, dotted, and dashed lines correspond to contributions of + , 1~ , and 2 + breakup, 
respectively. The experimental data are taken from Ref. 162 p. 



3.2. Differential breakup cross section for He scattering 

In Fig. U the calculated breakup cross sections are compared with the experimen- 
tal datdP on 6 He+ 12 C and 6 He+ 208 Pb reactions at 240 MeV/A. In the calculation, 
we take the same potentials as in Ref. [63D for n- 208 Pb and a- 208 Pb subsystems. The 
optical potential for n- 12 C subsystem is taken from the global nucleon-nucleus po- 
tential,^* while the optical potential for a- 12 C subsystem is constructed from the 
12 C+ 12 C potential at 200 MeV/AE3 by changing the radius parameter from 12 C to a. 
Nuclear breakup is dominant for 6 He+ 12 C scattering at 240 MeV/A, while Coulomb 
breakup to l - continuum is dominant for 6 He+ 208 Pb scattering. For 12 C target, 
the present theoretical result is consistent with the experimental data except for the 
peak of the 2 + -resonance around e = 1 MeV. Similar overestimations are also seen in 
the results of four-body distorted-wave Born approximation (DWBA). 66 ' For 208 Pb 
target, the present method underestimates the experimental data at e > 2 MeV. 
A possible origin of this underestimation is that the inelastic breakup reactions are 
not included in the present calculation. It has been reported in Ref. 166]) that the 
inelastic breakup effect is not negligible and thereby the elastic breakup cross section 
calculated with four-body DWBA underestimates the data. 

3.3. Summary on four-body CDCC 

We have proposed a new version of CDCC for treating four-body breakup. This 
method is called four-body CDCC. In the method, the three-body continuum of 
projectile is treated with the complex-scaling method. The validity of this method 
is checked for not only the elastic scattering but also the breakup reactions of 6 He. 
Clear convergence with respect to expanding the modelspace is seen in both the 
elastic and the breakup cross section. Four-body CDCC with the complex-scaling 
method is indispensable to study properties of unstable nuclei with two-neutron halo 
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structure. 

§4. Eikonal reaction theory 

In this section, we present an accurate method of treating neutron removal 
reactions at intermediate incident energies as an extension of CDCC and the Glauber 
model. This method is referred to as the eikonal reaction theory. This section is a 
brief review of Refs. |67|) . [68|) . 

4.1. Formulation 

Let us assume that a projectile (P) consists of a core nucleus (c) and a neutron 
(n). The scattering of P on a target (T) is then described by the three-body (c+n+T) 
Schrodinger equation 



~V% + h+U(r c ,r a )-E 



& = (4-1) 



with the interaction 



U = U™(r n ) + U^(r c ) + uP(r c ), (4-2) 



where h = T r + V(r) is the projectile Hamiltonian, R = (b, Z) stands for the 
coordinate between P and T, r x (x=n or c) represents the coordinate between x and 
A, and l4 N) and uP are the nuclear and Coulomb parts of the optical potential 
between x and T, respectively. Solving (|4-ip with the eikonal approximation, one 
can get the S-matrix operator 



S = exp 

for the operator 



iV \ dZO^UO 



(4-3) 



O = ~^e lKZ (4-4) 



with the wave-number operator K = ^j2[i{E — h) /h and the velocity operator 
i) = HK I fx of the relative motion between P and T, where V is the path order- 
ing operator. In the Glauber model, the adiabatic approximation is made as the 
secondary approximation in which h is replaced by the ground-state energy eo, and 
hence O^UO and V in (j4-3f) are reduced to U/(hvo) and 1, respectively, where vq is 
the velocity of P in the ground state relative to T. 

The operator O shows internal motions of c and n during the scattering. The 
movement is small for short-range nuclear interaction, but not for long-range Coulomb 
interaction. This means that uP is commutable with O with high accuracy: 

d^u^o ^UP/(hv ). (4-5) 

For the scattering of 31 Ne(lj>3/2) from a 208 Pb target at 240 MeV/nucleon, the error 
due to the replacement (|4-5[) is estimated with CDCC; it is 0.2% for the reaction 
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cross section <tr, 1.9% for the breakup cross section o"b u > 4.1% for the stripping cross 
section cr str . Using this replacement, one can get an important result 

S = S n S c (4-6) 

with 



S n = exp 
S c = exp 



^- I ~ dZU™ 



oo 
oo 



(4-7) 

/oo 
dZO\U^ +UP)6 . (4-8) 
-oo - 1 



Thus S can be separated into the neutron part S n and the core part S c . One can 
not calculate S c directly with Eq. (|4-8|) . because it includes the operators O and V. 
However, S c is the solution to the Schrodinger equation 



&c = 0, (4-9) 



■^Vl + h + uP(r c ) + uP(r c )-E 

when the eikonal approximation is made. One can obtain S c by solving Eq. (|4-9j) 
with eikonal-CDCO^ in which the eikonal approximation is made in the framework 
of CDCC. Non-eikonal corrections to S c can be easily made by solving Eq. (|4-9p 
with CDCC instead of eikonal-CDCC, although it is not necessary for the present 
intermediate scattering. As mentioned above, S n is obtained from Eq. (|4-7p . We 
can derive several kinds of cross sections with the product form ()4-6p . following the 
formulation on the cross sections in the Glauber modelP 1 '® 

4.2. One-neutron removal cross section for 31 Ne scattering 

The eikonal reaction theory is applied to one-neutron removal reactions for 
31 Ne+ 12 C scattering at 230 MeV/nucleon and 31 Ne+ 208 Pb scattering at 234 MeV/nucleon. 
Table [J presents several kinds of cross sections and the spectroscopic factor S = 
c^jf/c 1 ?!!- Thus 5[lp3/2] little depends on the target and less than 1, but 5[0/7/2] 
does not satisfy these conditions. In Ref. 131 j> . the Coulomb component of cj_ n [ 208 Pb] 
for a 208 Pb target is estimated to be 540 mb from the experimental values of 
<7_ n [ 208 Pb] and cr_ n [ 12 C]. In the eikonal reaction theory, the Coulomb component 
of (7_ n [ 208 Pb] agrees with o"b u [ 208 Pb] with good accuracy. The spectroscopic factor 
evaluated from the Coulomb component is 5' = 540/cj^ = 0.675 for the lp3/2 or- 
bit and 7.36 for the 0/7/2 orbit. Thus S' is consistent with S only for the lp3/2 
orbit. Hence, we can infer that the major component of the 31 Ne g . s . wave function 
is 30 Ne(0 + ) <g) lp3/2 (5 ~ 0.69). We adopt in the following this configuration. 

The potential V between c and n is not well known. Hence, S has a theoretical 
error coming from the potential ambiguity. The error is often estimated by changing 
the potential parameters by 30%. When the one- neutron separation energy B n of 
31 Ne is 0.33 MeV, S = 0.693 ± 0.133 ± 0.061 for 12 C target and 0.682 ± 0.133 ± 
0.062 for 208 Pb target, where the second and third numbers following the mean 
value stand for the theoretical and experimental uncertainties, respectively. Thus S 
includes a sizable theoretical error. This situation completely changes if we look at 
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Table I. Several kinds of cross sections and the spectroscopic factors for Ne-I- 1 C scattering at 
230 MeV/nucleon and 31 Ne+ 208 Pb scattering at 234 MeV/nucleon. The cross sections are 
presented in units of mb and the data are taken from Ref. I31j> . 







rz C target 






2m Pb target 






P3/2 


/r/a 


Exp. 


P3/2 


h/2 


Exp. 


OR 


1572.5 


1489.9 




5518.0 


4589.5 




C"bu 


23.3 


3.3 




799.5 


73.0 


(540) 


<m(-n) 


1463.5 


1458.6 




5151.5 


4524.2 




Cbu(-n) 


4.5 


1.0 




677.2 


60.5 




ffstr 


90 


29 




244 


53 






114 


32 


79 


1044 


126 


712 


S 


0.693 


2.47 




0.682 


5.65 





the asymptotic normalization coefficient (ANC) Canc- 69 ' When B Q = 0.33 MeV, 
Canc = 0.320 ±0.010 ±0.028 fin" 1 / 2 f or 12 C target and 0.318 ±0.008 ±0.029 fm" 1 / 2 
for 208 Pb target. Thus, Canc has much smaller theoretical errors than S. This 
means that the one-nucleon removal reaction is quite peripheral. 

4.3. The eikonal reaction thoery for two-neutron removal 

The eikonal reaction theory is applicable for two- neutron removal reactions. Here 
we show the application for the scattering of 6 He from 12 C and 208 Pb targets at 240 
MeV/nucleon. In this case, the projectile is a three-body system and hence four- 
body CDCC should be used. The potentials for n-target and a-target subsystems 
are calculated by the folding procedure shown in Sec. 12.21 The wave functions are 
obtained by spherical Gogny-HF calculations. Table [II] shows the integrated cross 
sections for two-neutron removal of 6 He. Our results are almost consistent with the 
experimental data!^ 1 Thus we can clearly see the reliability of the eikonal reaction 
theory for two-neutron removal reactions on both light and heavy targets. 

Table II. Integrated cross sections for two-neutron removal of 6 He. The cross sections are presented 
in units of mb and the experimental data are taken from Ref. I62[) . 





12 


C target 


208 


Pb target 




Calc. 


Exp. 


Calc. 


Exp. 


&ln str 


153.4 


127 ± 14 


353.6 


320 ± 90 


a 2n str 


29.0 


33 ± 23 


148.9 


180 ± 100 


C-2n 


198.5 


190 ± 18 


1016.6 


1150 ± 90 



4.4. Summary on the eikonal reaction theory 

We have presented an accurate method of treating neutron removal reactions at 
intermediate energies. In the theory, the nuclear and Coulomb breakup processes 
are accurately and consistently treated by CDCC without making the adiabatic 
approximation to the latter, so that the removal cross section calculated with the 
method never diverges even in the presence of the Coulomb interaction. For lower 
incident energies where the eikonal approximation is not perfectly accurate, one 
should make non-eikonal corrections to inclusive cross sections. This can be done 
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easily by using CDCC instead of eikonal-CDCC. 

Canc and S of the last neutron in 31 Ne are evaluated from the measured one- 
neutron removal reaction. Canc has a smaller theoretical error and weaker target- 
dependence than S. Thus, Canc is determined more accurately than S. When the 
last neutron of 31 Ne is in the lp3/2 orbit, S < 1 for B n < 0.6 MeV, and S and 
Canc have weaker target dependence. When the last neutron is in the 1/7/2 orbit, 
meanwhile, S > 1 and S and Canc have stronger target dependence. These results 
indicate that the last neutron is mainly in the lp3/2 orbit. This means that 31 Ne is 
deformed. This is consistent with our result in Sec. 12.21 The accuracy of the Glauber 
model is systematically investigated for deuteron scattering at 200 MeV/nucleon; see 
Ref. 68) for the details. 
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